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Abstract 

The role of image charges in nanoporous semiconductor materials is investigated within the 
framework of the effective mass and envelope function approximations. We show that nanomet- 
ric air bubbles in these materials can act as electron-trapping centers. This trapping capability 
originates from a deep stabilizing self-polarization potential well induced by the air - semiconduc- 
tor dielectric mismatch which can surpass the electroaffinity barrier. The trapping strength is a 
function of the pore size and the bulk parameters of the matrix material. A trapping parameter 
characteristic for each semiconductor material is defined. This parameter provides a simple way to 
ascertain the maximum pore size in a given material which is able to induce self-trapping of excess 
electrons. 

PACS numbers: 73.20.-r, 73.22.-f, 78.55.Mb 
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The progress made in silicon technology and the increasing miniaturization connected 
with a simultaneously increasing integration density require the development and integra- 
tion of a new generation of low dielectric constant materials.i^jMi^ The most promising 
strategy to diminish the effective dielectric constant of known materials relies on the 
introduction of pores in the material. This fact has raised vital interest in understanding 
the physical properties of porous materials . ^i^i^i^i^ Integration into microelectronic devices 
demands low pore interconnectivity and an extremely small pore size, thus forcing pores to 
be within the nanometric scale.— Pore sizes ranging from one to several tens of nanometers 
are nowadays commonly obtained. -^^^^ 

Due to the different dielectric response of the involved media, non-negligible surface 
image charges of free electrons (or holes) may arise at the interface between the air pore 
and the matrix material. These charges give rise to self-polarization of the inducing 
carriers and polarization of the Coulomb interaction between them.i^ The shape of the 
self-polarization potential is particularly sharp when the dielectric constant of one of 
the attached media is close to unity, yielding a narrow, deep stabilizing well near the 
interface that may eventually induce self-trapping of carriers . Both experimental and 
theoretical evidence of carrier trapping in dielectrically mismatched surfaces can be found 
in the literature . ^^I'^^i'^^i'^^i^^i^''' Since the largest effects of dielectric confinement take place 
in zero-dimensional nanostructures (quantum dots),— many works have been devoted to 
studying these systems under high dielectric mismatch conditions , ^i-^^i^^'^^i^^'^^i^'^'^^i^'-'i^^i^^i^^ 
particularly when the surrounding medium is air or a vacuum. However, to our knowledge, 
there are no studies that focus on the dielectric effects when the situation is reversed, that 
is, when the confined medium is air or a vacuum. 

In this paper we study the trapping capability of conduction band electrons in spherical 
nanopores (air bubbles) generated in different semiconductor materials. We find that, 
despite the barrier-acting nature of air cavities, dielectric confinement may induce electron 
localization in a narrow, deep attractive well within the pore, close to its surface. Depending 
on the pore size and the physical properties of the semiconductor we can encounter either 
electron trapping or scattering by the air bubbles, which leads to opposite implications in 
transport through nanoporous materials."^^ 
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The present study is carried out within the framework of the effective mass and envelope 
function approximations. Consequently, we employ a macroscopic treatment of Coulombic 
interactions, so that a parameter, the dielectric constant, characterizes the dielectric response 
of each involved material. The validity of such a treatment has been well-established for zero- 
dimensional semiconductor heterostructures similar to those presented here.— Our model 
consists of a spherical air bubble {Sair = 1, ^air — ^) surrounded by a semiconductor 
characterized by its dielectric constant e and conduction band electron effective mass m*. 
The corresponding one-particle effective mass Hamiltonian reads (in atomic units), 

where the first term is the generalized kinetic energy operator, and V^(r) represents the 
spatial confining potential. When the effective-mass approach and the envelope function 
approximation are used, the confining potential has a well defined step-like character at 
interfaces separating two different media, the rectangular steps being determined by the 
corresponding band offsets. In our case, as the mismatch involves air and a semiconductor, 
the band offset has been assumed to be equal to the semiconductor electroaffinity. Finally, 
05 (r) stands for the electron self-polarization potential coming from the image charges 
generated by the dielectric mismatch at the bubble edge. The calculation of this potential 
is carried out employing a dielectric function profile that changes smoothly within a thin 
interface -of the order of a lattice constant- between the air and the bulk semiconductor . 
It bypasses the (unphysical) self-polarization potential diverg^ces that arise at the interface 
when a step-like dielectric profile is employed (see Refs. yyilJllGllB for details). 
This model finds justification from the physical point of view, since the interface between 
two semiconductors (or between semiconductor and vacuum) is never perfectly sharp as 
the step-like model assumes. Conduction band ground state energies and wave functions 
were obtained by carrying out exact (numerical) integration of Eq. ([Q) by means of a 
discretization scheme on a grid extended far beyond the bubble radius. The explicit 
expression of 4> Jr) a nd a detailed description of the integration method employed can be 



found in Refs. 
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The bottom of Figure 1 shows a concrete example of the self-polarization potential 
profile for a -R=5 nm air bubble generated in a £=10 semiconductor material. As previously 
mentioned, the self-polarization potential yields a narrow, deep well within the nanopore, 
close to the interface, and has a weak repulsive character outside (see Fig. 1). The 
stabilizing effect of the self-energy well competes, then, with the destabilizing effect of 
the barrier that V^(r) presents inside the nanopore, thus yielding two different situations. 
On the one hand, if image charge effects dominate, the electronic ground state energy 
arises below the conduction band edge of the semiconductor, and its wave function is 
mainly localized within the nanopore (see Fig. 1). We will refer to these states as trapped 
states unless an amount of the total electron density higher than 30% spreads outside the 
nanopore. On the other hand, if V^(r) dominates, i.e. the semiconductor electroaffinity is 
large enough, the nanopore acts as a barrier, and the electronic density spreads throughout 
the semiconductor matrix. 

Figure 2 represents the maximum semiconductor electroaffinity Vmax that yields trapping 
of electrons vs. the natural logarithm of its dielectric constant. This is shown in Fig. 2a 
for a fixed value of the effective mass (m* = 0.5) and different values of the pore radius 
R ranging from 2 to 20 nm. It is also represented in Fig. 2b for a fixed value of the 
pore radius [R = 5 nm) and different effective masses of the semiconductor matrix rang- 
ing from 0.03 up to 1, this effective mass range accounting for most semiconductor materials. 

Several factors determine Vmax- On the one hand (see Fig. 2), as the semicon- 
ductor permittivity increases, Vmax also increases because the stabilizing well of the 
self-polarization potential becomes deeper, and thus raises the maximum value of the 
barrier height that admits electron trapping.^- The bubble radius also influences Vmax 
(see Fig. 2a) because the self-polarization potential scales as 1/R. Therefore, the 
larger the pore radius becomes, the weaker the trapping strength will be. Finally, the 
semiconductor effective mass m* also plays a role (see Figure 2b) because a part of the 
electronic density of a trapped electron leaks outside the bubble. Thus, a decrease in m* 
becomes an increase of the kinetic energy, i.e. an electronic energy destabilization. In other 
words, as m* diminishes, the trapping character of the air bubble is attenuated (see Fig. 2b). 
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An extra factor that can modify the trapping strength of a pore is a previously trapped 
electron. In order to analyze the magnitude of the Coulomb interaction (including polar- 
ization terms) and its influence on Vmax we have carried out CI j;^lculations on a system 
of two trapped electrons employing the same scheme as in Ref. The results obtained 
show that the total Coulomb energy contribution amounts to only about 0.1 eV, and, in 
consequence, the Vmax values are very slightly modified by the presence of a previously 
trapped electron in the air bubble. In other words, the associated Coulomb blockade effects 
exert a minor influence on the trapping strength of the air bubbles in comparison to the 
other variables R, e, m* and V . 

The lines shown in Fig. 2 represent linear regressions of the evaluated points with cor- 
relation coefficients above 0.99 in all cases. This linear behavior suggests the following 
exponential law relating Vmax vs. e, 

where Vq is the slope of the lines in Figure 2 and, as we show below, only depends on 
i?, thus accounting for the dependence of the trapping strength on the pore size. A is a 
second parameter (i?-independent) that accounts for the influence of the semiconductor 
effective mass on trapping. A comparison of the slopes of the (parallel) lines in Fig. 2b 
evidences the fact that Vq is m*-independent, while the change of slope vs. R for a fixed 
m* (Fig. 2a) reveals its dependence on R. The common x-intercept of lines in Fig. 2a 
reveals that A is i?-independent, while Fig. 2b shows the dependence of A on m*. Several 
other representations like Fig. 2a for different m* and Fig. 2b for different R (not shown) 
confirm the meaning and functional dependence of Vq and A. 

Vo has a simple decreasing linear dependence on i? in the range of bubble radii studied 
(see Fig. 3) while the m*-dependence of A that we find in the range of masses studied (0.03 
- 1) fits the following equation, 

A{m*) = 0.83 (m*)°-il (3) 
We can define a parameter t that is characteristic of a given semiconductor bulk and 
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which depends on the effective mass m*, permittivity e and elect roaffinitiy V as follows: 



t = (4) 

where A = A{m*) incorporates the m*-dependence of t. Equation (j3)) allows us to rewrite 
Eq. as, 

= *■ 

The trapping capability of a pore of radius R in a semiconductor matrix with a 
given m* and e depends on the electroaffinity V. The condition V = Vmax defines the 
border between trapping and scattering. At this limit, i.e. subject to the condition 
V = Vmax, t determines Vq and therefore the radius Rmax corresponding to this border. In 
other words, the trapping parameter t (an intrinsic property of the semiconductor bulk) 
determines the largest radius Rmax that a pore can reach without losing its trapping ca- 
pabilities. If the pore radius R is larger than Rmax it will act as a scattering center instead.— 



Figure 4 presents the relation between this trapping parameter t and Rmax, the maxi- 
mum pore radius allowing trapping capabilities. Several specific examples of well-known 
semiconductor materials, namely Si, Si02, and Ti02, are enclosed in the figure for the 
sake of illustration. The semiconductor parameters employed and the t values obtained 
are shown in Table I. We observe that, due to the large electroaffinity of Si, pores in this 
material are unable to localize the electronic density, regardless of their size. However, 
porous Si02 admits electron trapping in pores with radii up to about 13 nm. 



An especially interesting case is that of Ti02, due to the strong influence that temper- 
ature exerts on its permittivity.^^ Dielectric constants ranging from 6 to beyond 150 have 
been reported in the literature . ^^^^^ We consider e = 100 and e = 30 as being representative 
of two different temperature regimes (low and high, respectively). The results obtained 
(Fig. 4) show that at a low temperature {e = 100) electron trapping holds for bubble 
radii up to 7 nm. This trapping capability vanishes as the temperature increases (see 
the case e = 30 in Fig. 4). We emphasize that this dependence on temperature of Ti02 
as a trapping/scattering center could induce relevant electronic transport changes vs. 
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temperature in nanoporous samples of this semiconductor.— 

Summing up, we have investigated the possible behavior of nanopores (air bubbles) in 
semiconductor materials as carrier trapping centers. We found that the dielectric con- 
finement can be strong enough to overcome the barrier-acting nature of air, allowing the 
localization of electronic density within the nanopore. Large electron effective masses and 
dielectric constants, together with small electroaffinities and pore sizes, constitute the best 
scenario for electron trapping. A trapping parameter t accounting for the contribution that 
the specific semiconductor matrix makes to the trapping capability of its nanopores has been 
defined. This trapping parameter t provides a straightforward way to ascertain whether a 
specific nanoporous semiconductor material is able to induce electron trapping and, when 
this is the case, to obtain the largest radius of a trapping pore. 
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The employed effective mass m*=l for Ti02 may be found in Ref. 



effective masses, e.g. m*=10 (Ref. 



30|) or m*=20 (Ref. 



29l . However, much larger 



3lh have been also reported. A larger m* 



would yield a larger value for the parameter t, i.e., a larger trapping capacity. 
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Table 1. Semiconductor parameters employed and t values calculated (Eq. (jH)) for Si, Si02 
and Ti02. Two different dielectric constants are used to represent the regimes of low (T|) 
and high (T|) temperatures in Ti02. 





V{eV) m* 


e 


t{eV-') 


Si 


4.0 " 


0.26 ^ 


12 ° 


0.53 


Si02 


0.9 


0.5 ^ 




1.22 


Ti02 


3.9 


1.0 


100 (Ti) 


1.13 








30 (Tt) 


0.82 


" Ref. 


[22/ 


Ref. [22], 


"Ref. [2^ 


, Ref. [22] 
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FIG. 1: Top: ground state radial probability density of a trapped electron in an R=5 nm air- 
filled spherical cavity generated in an m*=0.5, e=10 and 1.8 eV electron-affinity semiconductor 
material. Bottom: self-polarization potential profile {(psir)) corresponding to this system. The 
vertical dotted line indicates the air bubble edge. 



FIG. 2: (a) Maximum semiconductor electroaffinity that allows trapping in nanopores as a func- 
tion of the natural logarithm of the semiconductor dielectric constant, calculated for m*=0.5 and 
different bubble radii. The lines displayed represent linear regressions of the calculated data. The 
evaluated points (19 in each series) were obtained with a precision of 0.1 eV, and have been omitted 
for the sake of clarity, (b) Same as (a) but for R =5 nm and different values of the effective mass 
of the semiconductor matrix. 



FIG. 3: Vo (see Eq. 0) vs. the air bubble radius R. The solid line represents the linear regression 
of the points displayed, with a correlation coefficient of 0.9974. The data employed corresponds 
to the case of m*=0.5, although we essentially obtain the same fit when other m* values are used 
(not shown). 



FIG. 4: Maximum pore radius Rmax that allows electron trapping as a function of the trapping 
parameter t of the semiconductor matrix. Example values for Si, Si02 and Ti02 have also been 
enclosed. 
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